Zero-Phonon Line Broadening and Satellite Peaks in Nanowire Quantum Dots: The 

Role of Piezoelectric Coupling 
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We investigate the influence of the one-dimensional character of the phonon modes in a catalyt- 
ically grown GaAs nanowire on the absorption spectrum of an embedded quantum dot, focusing 
on the contribution from the piezoelectric coupling. In general, the reduced dimensionality of the 
phonons leads to spectral side peaks and a zero-phonon line broadening due to the energetically 
lowest (acoustic) phonon mode. While the deformation potential predominantly couples to radial 
modes, the piezoelectric interaction can also couple strongly to modes of axial character, leading to 
additional absorption features. The contribution of the piezoelectric coupling to the zero-phonon 
line is negligible. 
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INTRODUCTION 

Recent epitaxial progress in the fabrication of catalyti- 
cally grown nanowires as well as the possibility to include 
heterostructures has led to an increased interest in quan- 
tum dots (QDs) in nanowires [l], 0]. QDs are interest- 
ing for the study of quantum information and cryptog- 
raphy d, 0, Q , requiring a good control of the dephasing 
mechanisms. Since the QDs are embedded in a strongly 
confined solid state matrix, it is important to take into 
account the influence of the geometry of the system on 
the phonon modes which offer an important decoherence 
channel. 

Recently, the influence of the phonons interacting with 
an electronic two-level system via the deformation po- 
tential coupling was investigated in a GaAs wire [6|. 
Here, we extend this analysis to include the interac- 
tion with phonons via the piezoelectric coupling which 
may be of particular relevance for wurtzite nanowires. 
In comparison with bulk modes which lead to continu- 
ous phonon side bands attached to the zero-phonon line 
(ZPL) 0, H, OH, the interaction with wire phonons 
leads to a set of phonon side peaks. These reflect the 
multitude of phonon modes arising due to the spatial con- 
finement in the lateral direction. In addition, while bulk 
phonons do not lead to a finite linewidth of the excitonic 
transition without resorting to higher order processes 
[1, EI El El, El, the lowest (acoustic) wire phonon mode 
leads to a temperature dependent ZPL broadening. In 
this paper, we focus on the piezoelectric interaction which 
can lead to additional absorption features as compared 
to the case of the deformation potential coupling due to 
different selection rules in the coupling mechanisms. 



following Ref. Q as well as the QD parameters used. 
The different electron-phonon coupling mechanisms are 
considered in Section . 

Hamiltonian: We consider an interband two-level sys- 
tem (conduction level |c), valence level \v)) interacting 
with the wire phonon modes via a diagonal coupling af- 
ter a 5-pulse excitation by a classical light field. After the 
initial excitation, the system is described by the Hamil- 
tonian 



H — H o + Hel- 



pii ! 



the first term, 



H = {e v a\a v + e c a\a c ) + hu qn b\ K b qK , 



(1) 



(2) 



describes the kinetics of the electronic two-level system 
and the phonons. Here, aj(^) denotes the creation (an- 
nihilation) operator of an electron in state \i) with energy 
£i, while b^ K (b qK ) denotes the creation (annihilation) op- 
erator of a phonon with quasi-momentum q and disper- 
sion u! qK , where the modes are labeled by k. The inter- 
action with the phonons is given by [H| 



cl— ph 



(q qK - q qK ) (b + fe f 



(3) 



Here, gf = J d 3 x v5*(x)V^ K (x)e i9 V 4 (x) is the electron- 
phonon coupling clement with the interaction potential 
^'( x ) = E 9K Vj K (x)e^(6 9K + 6% J, and the interaction 
is taken to be diagonal. No higher order phonon terms 
are considered here. 



THEORY 

In the following, we give a brief outline of the calcula- 
tion of the absorption spectrum and the phonon modes 



Absorption: The absorption is determined by the imag- 
inary part of the susceptibility via the macroscopic po- 
larization P(t) ~ d vc p(t) + c.c, where d vc is the in- 
terband dipole moment. The microscopic polarization 
p(t) = (aja c ) is calculated via the independent Boson 
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model (IBM) which can be solved analytically for an ar- 
bitrary number of phonon modes. The solution for t > 
is given by [IS El El 



p(t) = p(0)e- i(Wgap+A)t x 
J(lj 



(4) 



x exp 



duj- 



4n(uj) sin 2 y + 1 - e" iwt 



with the transition energy w gap = e c — e„, the Bosc- 
Einstein distribution = [exp^cj/fceT 1 ) — (tem- 
perature T) , and the polaron shift and the phonon spec- 
tral density, respectively given by 



A = 



dw 



\cJvv 



qn 12 



9K 



(5) 
(0) 



Phonon modes: The phonon modes are calculated 
within an isotropic continuum model based on Ref. [l6j |. 
Assuming free-surface boundary conditions, the solutions 
are uniquely determined, taking the general quantized 
form 



M,(x) 



1 



u qK (^ z (b qK +bl qK ) 



(7) 



(i = r,(p,z; N: number of unit cells). In the following 
discussions, we restrict to compressional modes (u v = 0). 
Torsional modes (u v ^ 0) do not couple via the de- 
formation potential coupling; for the piezoelectric cou- 
pling, this is taken to be an approximation. Furthermore, 
modes with a ^-dependence (so-called flcxural modes) 
are disregarded, since we assume a QD model which is az- 
imuthally symmetric, in which case the coupling to these 
modes vanishes. 

We consider a zinc-blende lattice in the (111) growth 
direction. The resulting lattice corresponds, in a 
nearest-neighbor approximation, to the wurtzite lattice 
171 ]. Calculations of strain in nanowires have shown 
that such a transformation is a good approximation 
[l8j |. The orientation is taken into account in the 
phonon modes via the velocities of sound which are now 
taken along the (111) direction. The dispersion rela- 
tion as well as two exemplary modes are shown in Fig. [1] 

Quantum dot model and wire parameters: For the 

QD model, we take a spherically symmetric charge dis- 
tribution with harmonic potential ^i(x) oc e~( r +z ^ 2tti , 
with Gaussian radii a c — 5.8 nm and a v = 3.19 nm. For 
the wire, we use GaAs parameters: deformation poten- 
tials D c = -14.6 eV, D v = -4.8 eV mass density p m 
— 5370 kg/m 3 , longitudinal and transverse sound veloc- 
ities vi — 5400 m/s, Vt — 2800 m/s, dielectric constant 
£ s = 12.9, and wire radius R — 25 nm. The piezoelectric 
constants for wurtzite GaAs are not known. We thus cal- 
culate the constants from the zinc-blende value ei4, again 
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FIG. 1: (left) Dispersion relation and (right) two exem- 
plary modes for q = 0.002 nm -1 of the one-dimensional wire 
phonons. 



using the transformation described above [17j • Using the 
value ei4 = —0.16 Cm~ 2 [lj|, we obtain an estimate 
for the wurtzite piezoelectric constants eis — e^i = 0.09 
Cm -2 and 633 = —0.18 Cm -2 . We note that larger values 
for ei4 can be found in the literature [2(| and that the 
transformation formula is an approximation, and thus 
the effect of the piezoelectric coupling could certainly be 
enhanced compared to the results shown here. 



ELECTRON-PHONON COUPLING 
MECHANISMS 

We consider two electron-phonon coupling mecha- 
nisms: the deformation potential and the piezoelectric 
coupling. In the following, we discuss and compare the 
coupling elements and the resulting phonon spectral 
densities. We restrict the discussions to the six energeti- 
cally lowest phonon modes. 

Deformation potential: The coupling of the electronic 
system to the phonons via the deformation potential is 
just given by the divergence of the displacement 



15] 



V 



Di V • u(x) 



(8) 



where Di is the deformation potential. The calculated 
coupling constants are shown in Fig. [Deleft), while the 
corresponding spectral density is shown in Fig. fright). 
Modes 3 and 6 are almost purely radial for q w (cf. 
Fig. Q]), thus they couple strongly via the deformation 
potential. However, due to the density of states, the 
spectral density of mode 3 is relatively weak. The axial 
modes 4 and 5 only couple weakly via the deformation 
potential for q rj due to the restriction to diagonal 
strain elements Sij — (diUj + djUi)/2 in Eq. ([8]). 

Piezoelectric coupling: For the piezoelectric coupling, 
the displacement strain causes an electrostatic potential 
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FIG. 2: (left) Coupling elements \g«% - g%\ 2 {q) (L = 1 nm) 
and (right) phonon spectral density J(ui) for the deformation 
potential coupling for the energetically lowest phonon modes. 



FIG. 4: Comparison of (left) the phonon spectral densities 
J(w) and (right) the absorption spectra at T = 77 K for the 
deformation potential and piezoelectric coupling. 



which is determined via 

£ £ s A^(x) = V • P pz , 



(9) 



where P p = ek,ijSij (piezoelectric constants e-k,ij) 
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22L |23j. For the wurtzite structure, only three in- 
dependent piezoelectric constants remain, leading to the 
polarization 



P pz (x) = eis (d r u z (x) + d z u r (x)) e r 



+ <j e 3 i ( <9 r w r (x) + -u r {x) ) + e 33 d z u z (x) 



(10) 



Assuming that the potential vanishes for r — > oo using 
e s = 1 outside the wire, we obtain an analytical solution 
for the electron-phonon potential V(x) = — e</>(x). 

In Figs. EJleft) and Slight), the coupling elements 
and the corresponding spectral density is shown, respec- 
tively. Compared to the deformation potential coupling, 
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ABSORPTION SPECTRA 

We now calculate the absorption within the IBM 
for the one-dimensional phonon modes. In Fig. [5J the 
absorption spectrum for (left) a QD interacting with 
the wire modes is compared to (right) a QD exposed to 
bulk phonons for different temperatures. In both cases, 
a phcnomcnological (radiative) dephasing of T 2 = 1 /ieV 
is used [7|, and the system is excited resonantly on the 
polaron shifted transition energy E = Ti{uj &av + A). In 
Fig. Upright), the different coupling mechanisms are 
compared in the absorption for the wire modes. 
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FIG. 5: Absorption spectrum of a QD interacting with (left) 
wire phonons via the deformation potential and the piezoelec- 
tric coupling and (right) bulk phonons via the deformation 
potential coupling for different temperatures. 



FIG. 3: (left) Coupling elements \g%« - g^\ 2 (q) (L = 1 nm) 
and (right) phonon spectral density J{uj) for the piezoelectric 
coupling for the energetically lowest phonon modes. 

modes 4 and 5 couple strongly here for q ss since the 
piezoelectric potential [Eq. (fTO)) ] takes into account the 
nondiagonal strain elements Sij which favor axial modes 
for small q (cf. Fig. [I}. This interplay between the de- 
formation potential and the piezoelectric coupling can be 
seen nicely when comparing the spectral densities of the 
two coupling mechanisms, shown in Fig. [5Jleft). 



Side peaks: Compared to the case of bulk phonons, 
where the linear dispersion of the acoustic phonon mode 
leads to continuous side bands attached to the ZPL, 
the discrete wire modes lead to a series of phonon 
side peaks. The side bands are strongly temperature 
dependent, where higher temperatures lead to both 
increased dephasing as well as a broader spectrum due 
to multi-phonon scattering processes. Furthermore, at 
low temperatures, the absorption/emission asymmetry 
can be seen. The peaks in the absorption spectra 
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correspond to the extrema in the dispersion, where the 
density of states has a singularity, whereas their strength 
is determined by the magnitude of the coupling at the 
extrema [f|. In addition to the radial modes 3 and 6 
which couple strongly via the deformation potential, the 
axial modes 4 and 5 can be seen in the absorption which 
interact via the piezoelectric coupling. The dependence 
of the scattering mechanism on the mode character can 
be seen more clearly in Fig. Upright), where the two 
mechanisms are compared. 

ZPL broadening: As was shown in Ref. [6|, the phonon 
mode extending to vanishing energy leads to a tempera- 
ture dependent broadening of the ZPL. This can be seen 
nicely in Fig. [Deleft) for increasing temperature. For 
the piezoelectric coupling, no considerable contribution 
to the ZPL broadening is found due to the form of the 
piezoelectric potential for q — > 0. 

CONCLUSION 

We have studied the interaction of a quantum dot 
in a catalytically grown GaAs nanowire via the defor- 
mation potential and piezoelectric electron-phonon cou- 
pling. In general, the electron-phonon interaction leads 
to discrete side peaks, which are strongly temperature de- 
pendent, and a broadening of the zero-phonon line. Since 
the deformation potential couples strongly only to radial 
modes, the piezoelectric coupling can become important 
for modes which are mainly axial and thus can play an 
important role for the determination of the optical spec- 
trum. The inclusion of the piezoelectric coupling does not 
lead to an appreciable contribution to the zero-phonon 
line broadening. It should be mentioned that larger val- 
ues of the piezoelectric constants than used here have 
been reported in the literature, and thus the influence 
of the piezoelectric coupling could become more impor- 
tant. Furthermore, the piezoelectric coupling is strongly 
dependent on the quantum dot geometry, and thus the 
relative as well as the absolute strengths of the coupling 
constants must be considered for each case. 

As an outlook, the size and shape dependence of 
the piezoelectric coupling to the one-dimensional wire 
phonon modes must be further investigated. In addition, 
the torsional modes, which are neglected in this investiga- 
tion, should be considered to ascertain their interaction 
strength via the piezoelectric coupling. 
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